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Abstract 

We study a new flavor symmetric model with non-Abelian discrete symmetry T13. The T13 
group is isomorphic to Z13 x Z3, and it is the minimal group having two complex triplets in 
the irreducible representations. We show that the X13 symmetry can derive lepton masses and 
mixings consistently. Moreover, if we assume a gauge-singlet fermionic decaying dark matter, its 
decay operators are also constrained by the T13 symmetry so that only dimension six operators 
of leptonic decay are allowed. We find that the cosmic-ray anomalies reported by PAMELA and 
Fermi-LAT are well explained by decaying dark matter controlled by the T13 flavor symmetry. 
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1 Introduction 



Despite the great success of the Standard Model (SM) of the elementary particle physics, the ori- 
gin of flavor structure, masses and mixings between generations, of matter particles are unknown 
yet. In order to overcome these problems, plenty of models based on the principle of symmetry, 
flavor symmetry, have been discussed. In particular, the fact that the lepton mixing matrix (Maki- 
Nakagawa-Sakata matrix) Umns shows very good agreement with the tri-bi maximal form [1] implies 
that flavor structure is originated from a symmetry. Among them, non-Abelian discrete symmetries 
are well discussed as plausible possibilities [2J. 

On the other hand, it has been established that about 23 % of energy density of the universe 
consists of Dark Matter (DM) [3]. Indirect detection experiments of DM, PAMELA [I] and Fermi- 
LAT [5j [6], reported excess of positron and the total flux (e + + e~) in the cosmic ray. These 
observations can be explained by scattering and/or decay of TeV-scale DM particles. Since PAMELA 
measured negative results for anti-proton excess [TJ, leptophilic DM is preferable. 

Even so, if the main final state of scattering or decay of DM is t + t~, this annihilation/decay 
mode is disfavored because it will overproduce gamma-rays as final state radiation [8]. This may 
indicate that if the cosmic-ray anomalies are induced by DM scattering or decay, these processes also 
reflect flavor structure of the theory. There are several papers in which the DM nature is related to 
flavor symmetry P [10l HH [121 H31 HH HSl US]. 

While there are many models for the DM, we consider decaying DM model in this paper [T7J [181 
[T9] . In such scenarios, no excess of anti-proton in the cosmic ray [7J implies that lifetime of the DM 
particle should be of O(10 26 ) sec. This long lifetime is achieved if the TeV-scale DM (gauge singlet 
fermion X) decays into leptons by dimension six operators LELX/K 2 suppressed by GUT scale 
A ~ 10 16 GeV [18J. In this case, the lifetime of the DM is estimated as T" 1 ~ ((TeV) 5 /A 4 ) -1 ~ 10 26 
sec. However, in general, there are gauge invariant dimension four decay operators which induce 
rapid DM decay, and several dimension six operators which induce DM decay into quarks, Higgs and 
gauge bosons. Therefore one has to solve at least two problems for decaying DM models: i) why the 
lifetime of the DM is so long? ii) why the DM decays mainly into leptons? In ref. [9], we and our 
collaborators have shown that A 4 flavor symmetry can allow only the LELX/K 2 operator and forbid 
the other undesirable operators, and determine flavor structure of DM decay mode. In that model, 
the symmetry allows only flavor universal leptonic decay mode, and the cosmic-ray anomalies are 
explained by fermionic DM decay. 

In this paper, we show that similar argument is possible in T 13 flavor symmetry model as an exten- 
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sion of the model. The Ti 3 group, isomorphic to Z13 x Z3 group, is non-Abelian discrete subgroup 
of SU(3). For the lepton sector, tri-bi maximal form can be derived by embedding three generations 
into triplet representations of flavor symetries. In this point of view, A4 is the minimal group which 
has a triplet, and the T 13 group has two complex triplets in the irreducible representations. However, 
since multiplication rules of the T 13 group is very different from those of A4, we have a new texture 
of mass matrices of the lepton sector. As a result, although the A4 model requires SU{2) L triplet 
Higgs bosons A with heavy mass ttt-a and small vacuum expectation values (VEVs) v a to make the 
leptonic mixings and to suppress additional dimension five DM decay operator HA^LX, our Ti 3 
model does not require A. The DM decay operators are also constrained by the T 13 symmetry so 
that only LELX/K 2 operators are allowed like v4 4 model of ref. [9]. However unlike the A 4 model, 
DM decay mode depends on mixing matrices in general. In this paper we choose a particular set of 
parameters of the lepton sector, and show that the cosmic-ray anomalies can be well-explained by 
fermionic DM decay controlled by T 13 symmetry. 

This paper is organized as follows. We briefly discuss group theory of T 13 and lists the multipli- 
cation rules in the next section. In the section 3, we construct mass matrices of the lepton sector in 
definite choice of Ti 3 assignment of the fields, and show that there exists a consistent set of param- 
eters. In the section 4, we show that only desirable dimension six DM decay operators are allowed 
by T 13 symmetry and that leptonic decay of the DM by those operators shows good agreement with 
the cosmic-ray anomaly experiments. The section 5 is devoted to the conclusions. 

2 T13 group theory 

First of all, we briefly review the non-Abelian discrete group Ti 3 , which is isomorphic to Z13 x Z 3 
[2"U| l2~Tj . The T 13 group is a subgroup of SU(3), and known as the minimal non-Abelian discrete 
group having two complex triplets as the irreducible representations. We denote the generators of 
Z13 and Z 3 by a and b, respectively. They satisfy 
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1, ab = ba . 



(2.1) 



Using them, all of 7\ 3 elements are written as 
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Table 1: Characters of Ti 3 . & is defined as the complex conjugate of £j. 



The generators, a and 6, are represented e.g. as 
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where p = e 2j7r//13 . These elements are classified into seven conjugacy classes, 
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(2.3) 



(2.4) 



The T 13 group has three singlets with = 0, 1, 2 and two complex triplets 3i and 32 as 
irreducible representations. The characters are shown in Table [TJ where £i = p + p 3 + p 9 , £2 = 
p 2 + p 5 + p 6 , and u = e 2in/3 . 

Next we show the multiplication rules of the T 13 group. We define the triplets as 
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where the subscripts denote Z13 charge of each element. 
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The tensor products between triplets are obtained as 
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The tensor products between singlets are obtained as 
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The tensor products between triplets and singlets are obtained as 
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In the following section, we discuss mass matrices of the lepton sector determined by the T 13 
flavor symmetry. 



3 Lepton masses and mixings 

In this section, we discuss the lepton masses and mixings in the setup shown in Table [2j Here, Q, U, 
D, L, E, H(H') and X denote left-handed quarks, right-handed up-type quarks, right-handed down- 
type quarks, left-handed leptons, right-handed charged leptons, Higgs bosons, and gauge singlet 
fermion, respectively. We introduce two T 13 triplet Higgs bosons H(3i) and H(3 2 ) which couple to 
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leptons, and three T i3 singlet Higgs bosons H'(l ^ 2 ) couple to quarks. This is realized by appropriate 
choice of an additional Z 3 symmetry. Since we concentrate on the lepton sector in this paper, T 13 
charge of quarks and H' are assigned to the singlets. Therefore, mass matrices in the quark sector 
are not constrained, while those in the lepton sector are determined by the T 13 symmetry. A gauge 
and T 13 singlet fermion X is also introduced in addition to the SM matter fermions. Since the 
Yukawa couplings LXH are forbidden by the T 13 symmetry, this new fermion X does not work as 
right-handed neutrino, and left-handed Majorana neutrino mass terms are generated by dimension 
five operators LHLH . The fermion X is a DM candidate decaying into leptons by dimension six 
operators LELX. However, we postpone the discussion on this point to the next section. 

For the matter content and the T 13 assignment given in Table [21 the charged- lepton and neutrino 
masses are generated from the 7\ 3 invariant operators 



C Y = V2a e ELH c {3 2 ) + V2b e ELH c {3 1 ) 



+ j-LH(3 2 )LH(3 2 ) + jLH(3 2 )LH(3 1 ) + j-LH{3 2 )LH{3 1 ) + h.c, 



(3-1) 



where H c = eH* and the fundamental scale A = 10 11 GeV. After the electroweak symmetry breaking, 
the Lagrangian Eq.f l3.ip gives rise to mass matrices of charged leptons M e and neutrinos M v 
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where the VEVs are defined as 

(3iy> = ^ (H(3 2 y) = E M + ^) = ( 246 GeV ) 2 - 

Now we give a numerical example. By the following choice of parameters 

V! = 0.118297 GeV, v 2 = 179.257 GeV, v 3 = 1.99413 GeV, 
«i = 10 GeV, v 2 = 168.164 GeV, v 3 = 0.0483626 GeV, 
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Table 2: The T13 and Z3 charge assignment of the SM fields and the dark matter X, where oj = e 2j7r / 3 . 
the mass matrices Eqs. (13. 2\i and (13.31) give rise to mass eigenvalues and related observatives as 



Am 2 21 
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/ 



which are all consistent with the present experimental data [22]. In particular in the case of U e i = 1, 
the mass matrices Eqs. (13. 2p and (13. 3p require normal hierarchy m v \ < m U 2 < m v3 of the neutrino 
masses and U^ NS 7^ 0. A comprehensive analysis of the T13 symmetry models will be published 
elsewhere [23] . In the following analysis, we assume these parameters to discuss decaying dark matter. 



4 Decaying dark matter in the T13 model 

The cosmic-ray anomalies measured by PAMELA [4] and Fermi-LAT [51 [6] can be explained by DM 
decay with lifetime of order ~ 10 26 sec. If the DM decays into leptons by dimension six operators 
LELX/A 2 , where X is the gauge singlet fermionic DM, such long lifetime can be achieved. In general, 
however, there exist several gauge invariant decay operators of dimension four which induce rapid 
DM decay, and of dimension six which include quarks, Higgs and gauge bosons in the final states. In 
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this section, we first show that the T13 flavor symmetry forbids those undesired operators. Next we 
show by calculating the positron (electron) flux that the scenario, which has a particular generation 
structure of DM decay vertices, is possible to excellently describe the cosmic-ray anomalies. 

4.1 Dark matter decay operators 



Dimensions 


DM decay operators 


4 


LH C X 


5 




6 


Lelx, WhLh c x, (h c Yd^h c eyx, 




QDLX, UQLX, LDQX, U^DE^X, 




D^^D^LX, D»D^H C LX, 




B^La^H c X, W^La^T a H c X 



Table 3: The decay operators of the gauge-singlet fermionic dark matter X up to dimension six. B^, W^ u , 
and are the field strength tensor of hypercharge gauge boson, weak gauge boson, and the electroweak 
covariant derivative. 

In our model, as mentioned before, a gauge-singlet fermion X is introduced as the dark matter 
particle in addition to the SM fermions. By assuming that the baryon number is preserved at least 
at perturbative level, it turns out that there exist various gauge invariant operators up to dimension 
six shown in Table [3] [24] . From the table, one finds that the dark matter X can in general decay 
into not only leptons but also quarks, Higgs, and gauge bosons at similar rates by dimension six 
operators. Furthermore, a rapid decay of DM is induced if the dimension four Yukawa operator 
LH C X is allowed. One may try to impose an Abelian (continuous or discrete) symmetry to prohibit 
unwanted decay operators, but it is shown in ref. [9] that it does not work but flavor symmetry 
does. However in the A± model discussed in ref. [9], SU(2)l triplet Higgs bosons A are introduced 
in order to give the mixings of the lepton sector, and these give rise to dimension five DM decay 
operators. In order to avoid rapid DM decay by the dimension five operators, small VEVs (A) and 
large mass wa are required. On the other hand, lepton masses and mixings are generated only by 
SU(2)l doublet Higgs bosons in the T 13 model as shown in the section 3. 

Remarkably, by the field assignment of Table |2j all the decay operators listed in Table |3] except 
for LELX are forbidden due to this single symmetry. Consequently the DM mainly decays into three 
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leptons. With the notation Li = {y^ti) = (U e L)ia(v a ,£ a ) and Ei = (U e R)ipEp (i = 1,2,3, a, (3 
e, /i,r), the four- Fermi decay interaction is explicitly written as 

A 3 



£decay = ^^(LiE^LjX + h.C. 
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i=l 
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As seen from Eq. (14. ip . decay mode of the DM particle X depends on the mixing matrices U e (L,R), 
which are given in Eq. (13. 7\i . 



4.2 Positron production from dark matter decay 

Next, we consider the branching fraction of the DM decay through the T 13 invariant Lagrangian 
Eq. (l4.ip . Due to the particular generation structure, the dark matter X decays into several tri- 
leptons final state with the mixing-dependent rate. The decay width of DM per each flavor (T Q ^ = 
F(X — > v a £p£~)) turns out to be 
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where mx is the DM mass, and 



U, 



a/3"/ 



(4.2) 



(4.3) 



Here we have omitted the masses of charged leptons in the final states. The flavor dependent factor 
U a /3^ gives a factor three if one takes the sum of flavor indices a, (3 and 7. Therefore, the branching 
fraction of each decay mode is given by 

BR(X uJ+£~) = i (U afh + U aif} ) . (4.4) 

The DM mass mx and the total decay width T = J2 a p 7 r aj g 7 are chosen to be free parameters in 
the following analysis. 

Given the decay width and the branching fractions, the positron (electron) production rate (per 
unit volume and unit time) at the position x of the halo associated with our galaxy is evaluated as 

~dN„± 



Q(E,x) = n x (x)Tj2MX ->• /) 

/ 



dE 



(4.5) 
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where [dN e ±/dE]f is the energetic distribution of positrons (electrons) from the decay of single DM 
with the final state '/'. We use the PYTHIA code [25] to evaluate the distribution [dN e ±/dE]f. The 
DM number density nx(x) is obtained by the profile p(x), the DM mass distribution in our galaxy, 
through the relation p(x) = rrixnx{x). In this work we adopt the Navarro- Frank- White profile [26] . 




where p Q ~ 0.30 GeV/cm 3 is the local halo density around the solar system, r is the distance from 
the galactic center whose special values r ~ 8.5 kpc and r c ~ 20 kpc are the distance to the solar 
system and the core radius of the profile, respectively. 

As seen from Eqs. (13.71) and f!4.2p . the DM decays into r as well as e and p ± in the equal rate. 
However, in the final states decay into hadrons, and such hadronic decays are suppressed by the 
electroweak coupling and the phase space factor. As a result, pure leptonic decays give dominant 
contributions, and it is consistent with no anti-proton excess of the PAMELA results [7]. On the 
other hand, the injections of high-energy positrons (electrons) in the halo give rise to gamma rays 
through the bremsstrahlung and inverse Compton scattering processes. These gamma-ray flux [27] 
from leptonically decaying DM may constrain mass and lifetime of the DM [281 122] • We will discuss 
this point later. Therefore, we concentrate on calculating the e ± fluxes in what follows. We follow 
ref . [9] for diffusion model describing the propagation of positrons and electrons [301 ED E2] , and 
backgrounds [301 E3] ■ 



The positron fraction and the total flux [$ e -] total + [$e+] total are depicted in Figure [T]for the scenario 
of the leptonically decaying DM with T 13 symmetry. For the DM mass mx = 2, 2.5, and 3 TeV, the 
results are shown with the experimental data of PAMELA and Fermi-LAT. The total decay width 
T is fixed for each value of DM mass so that the best fit value explains the experimental data. With 
a simple \ 2 analysis, we obtain T _1 = 8.4 x 10 25 , 6.95 x 10 25 , and 5.9 x 10 25 sec for m x = 2.0, 2.5, 
and 3.0 TeV, respectively. One can see that in the mx = 2.5 TeV case, both experiments are well 
explained in the Ti 3 model. Recent studies [281 129] suggest that mass and lifetime of the decaying 
DM are strongly constrained by gamma-ray measurement from cluster of galaxies, and that allowed 
region which can simultaneously explain the PAMELA and Fermi-LAT results does not exist. To 
avoid these constraints, mass and lifetime of the DM should be lighter than ~ O (TeV) and longer 
than ~ O(10 ) sec, respectively. In that case, only the PAMELA results can be explained by the 
decaying DM. 




Pnfw [x 



Pe 



r (r Q + r e ) 2 
r(r + r c ) 2 



(4.6) 



4.3 Results for PAMELA and Fermi-LAT 
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Figure 1: The positron fraction [3] and the total e + + e~ flux [5j |6] predicted in the leptonically-decaying 
DM scenario with T13 symmetry. The DM mass is fixed to 2.0, 2.5, and 3.0 TeV. As for the DM decay 
width used in the fit, see the text. 



5 Conclusions 

We have considered a new flavor symmetric model based on a non-Abelian discrete symmetry T13. 
The T13 group, isomorphic to Z13 x Z3, is the minimal group which contains two complex triplets in 
the irreducible representations. The form of mass matrices are determined by the assignment of T 13 
charges and the multiplication rules. We have shown that masses and mixings in the lepton sector 
are derived in the T13 model consistently. Thanks to the complexities of the T13 group compared to 
A4, both the leptonic masses and mixings are made only by SU(2)l doublet Higgs bosons. 

We have also shown that the decay of gauge-singlet fermionic dark matter can explain the cosmic- 
ray anomalies reported by the PAMELA and Fermi-LAT experiments. It is known that if the dark 
matter is TeV-scale fermionic particle, its longevity of order 10 26 sec can be derived from dimension 
six four-fermi operators suppressed by a large scale of new physics. The T13 symmetry forbids DM 
decay of final states with quarks, Higgs and gauge bosons, and allows only leptonic decay. Moreover, 
it determines the DM decay mode so that tauon final state does not give dominant contribution. 
We found that due to the fermionic DM decay controlled by the T13 flavor symmetry, the cosmic-ray 
anomalies are well-explained. 

In this paper, we have explicitly given a numerical example of one consistent set of parameters in 
the mass matrices of the lepton sector. For completeness, a comprehensive analysis of mass matrices 
and its phenomenology in T13 symmetric models will be published elsewhere |23j . 
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